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We theoretically consider a Josephson junction formed by a ferromagnetic spacer with a strong 
spin-orbit interaction or a magnetic spin valve, i.e., a bilayer with one static and one free layer. 
Electron spin transport facilitates a nonlinear dynamical coupling between the magnetic moment 
and charge current, which consists of normal and superfluid components. By phenomenologically 
adding reactive and dissipative interactions (guided by structural and Onsager symmetries), we 
construct magnetic torques and charge pumping, whose microscopic origins are also discussed. A 
stability analysis of our coupled nonlinear systems generates a rich phase diagram with fixed points, 
limit cycles, and quasiperiodic states. Our findings reduce to the known phase diagrams for current- 
biased nonmagnetic Josephson junctions, on the one hand, and spin-torque driven magnetic films, on 
the other, in the absence of coupling between the magnetic and superconducting order parameters. 



PACS numbers: 72.25.-b,74.50.+r,74.20.Rp,75.70.Cn 

Hybrid structures with ferromagnet (F)|normal-metal 
(N) interfaces have garnered much attention over the past 
few decades owing to their application in spintronic de- 
vices. Injecting a spin current into such a system ex- 
erts a torque on the magnet pQ, which can induce pre- 
cession and even reversal [2J, allowing for manipulation 
of the magnetic order parameter in nanoscale structures 
without an external magnetic field [3]. Because of the 
nonlinear nature of the ensuing magnetic dynamics, such 
devices offer observation of effects traditionally seen in 
nonlinear dynamical systems: Phase locking, hysteresis, 
bifurcations, and chaos are readily observed [3]. 

In consideration of a superconductor (S)|F|S het- 
erostructure, one may expect the Josephson effect to be 
suppressed due to the rapid decay of a singlet pair in- 
side the ferromagnet. Recent experiments [5], however, 
observed superconducting transport through a strong fer- 
romagnet between two conventional (s-wave) supercon- 
ductors. With the expectation that the triplet compo- 
nent of the superconducting condensate can penetrate 
long lengths into a ferromagnet, the preservation of this 
signal suggests a spin singlet-to-triplet conversion at the 
interfaces [6] . The unexpected persistence of a supercur- 
rent through the magnet forecasts a new kind of spin- 
tronic device that manipulates the Josephson junction 
by the ferromagnet and, conversely, ferromagnetic layer 
by the superconducting condensate [?_H§]. 

Previous analyses [71-110] have considered equilibrium 
interactions between magnetic and superconducting or- 
der parameters, which naturally induce a reactive cou- 
pling. In contrast, in our description, we introduce 
nonequilibrium interactions consistent with the symme- 
tries of the structure and obeying Onsager reciprocity 
[llj . This treatment allows the addition of both dissi- 
pative and reactive couplings between the magnet and 
superconductor that may in practice be crucial in the un- 
derstanding of ferromagnetic Josephson junctions, anal- 
ogous to the importance of Slonczewski [1 -3 and spin- 
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FIG. 1. Schematics of our magnetic Josephson junctions. 
The directions of junction layering, applied current J, internal 
Rashba field E (a), and direction of the static ferromagnetic 
layer m 3 (b) all lie along the z axis. <j> is the phase difference 
between the superconducting leads. 



pumping |12j terms in the theory of spin-transfer torques. 
Such effects cannot be fully captured by quasiequilibrium 
free-energy considerations. We expect the dissipation to 
be governed by the quasiparticle excitations in the super- 
conductors in concert with the microscopic processes in 
the ferromagnet (such as magnon-phonon and magnon- 
magnon interactions in insulators) that are responsible 
for their Gilbert damping (which, in turn, is known to 
persist down to very low temperatures). 

In order to provide specific examples, we consider (a) 
an S|F|S hctcrostructure with a Rashba spin-orbit inter- 
action (SOI) in a thin ferromagnetic interlayer (neglect- 
ing the vector potential and associated phase shift caused 
by its magnetic moment [5]) and (b) S|F|N|F|S het- 
erostructure wherein one ferromagnetic layer is pinned, 
m s , and the other free, m/. See Fig. [I] The corre- 
sponding spin-dependent Hamiltonians mix the singlet 
and triplet superconducting components [6] , allowing the 
superfluid to penetrate into the magnet and exert spin 
torque and carry spin pumping (since a triplet Cooper 
pair is a spin-1 object) that are analogous to those as- 
sociated with normal quasiparticles (spin-1/2 objects). 
In particular, as a simple model to demonstrate proof 
of concept, we take the device geometry to be rotation- 
ally symmetric along the axis associated with the Rashba 
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interaction, as sketched in Fig.jlja), and along the direc- 
tion of the fixed magnetic layer, as sketched in Fig. |ljb) . 
By analyzing the stability and dynamics of our model, 
we outline a phase diagram of the coupled system as a 
function of applied magnetic field and current bias. 

The phenomenological equation of motion of an iso- 
lated ferromagnet sufficiently well below the Curie tem- 
perature is given by the Landau-Lifshitz-Gilbert (LLG) 
equation [T3] 

m = —7m x H + am x m , (1) 

where 7 is the gyromagnetic ratio and a is the dimension- 
less Gilbert damping. We use a normalized form of this 
equation, in which the (unit) magnetic direction vector 
m = M/M,, M s = |M| (saturation magnetization), is 
dimensionless. H = — V~ 1 dF/dM is the effective mag- 
netic field and F, M, and V are the free energy, magneti- 
zation vector, and volume, respectively. In the spin-valve 
model, Fig. [TJb), m will denote the free layer m/. 

We consider the resistively-shunted junction (RSJ) 
model for the Josephson junction, wherein the device 
is composed of conventional superconductors with some 
Ohmic conductance a in the junction [14] . Additionally, 
we take the capacitance to be zero, which precludes RC- 
type delays in the coupled dynamics. The corresponding 
Josephson relations (for a static magnetization) are 

Q = I c sin 4> + (jV , j>= (2) 

where V is the voltage drop across the junction. (f> is the 
phase difference between the superconducting reservoirs 
and Q is the charge transported by the junction. The 
supercurrent is proportional to the critical current, I c — 
(2e/h)Ej, where Ej parametrizes the Josephson energy 
—Ej cos 4>. We note that Eq. ([2| is dictated by gauge 
symmetry and, in anticipation of the arguments to follow, 
is a manifestation of Onsager reciprocity in the dynamics 
of Q and (j). 

Under time reversal, m — > m and am x m — > —am x 
m. The term proportional to a in the LLG equation 
thus reflects irreversible processes. We characterize such 
terms as dissipative. <r, likewise, parametrizes Ohmic 
dissipation of normal fluid. All other terms thus con- 
sidered so far are reactive. Couplings between the free 
ferromagnet and superconductor at the level of the free 
energy, induced by the static magnetic layer or SOI, are 
restricted by the symmetries of our structure. Our device 
geometries shown in Fig. [T] are invariant under rotations 
about the z axis as well as certain combinations of dis- 
crete symmetries. Because both the exchange interaction 
between the magnetic layers of our spin- valve device and 
the Josephson energy are individually preserved under 
the symmetries of the combined system, the product of 
these interactions must also be permitted [7J. However, 
the interlayer F|N|F spin- valve exchange is usually very 



small (except for the thinnest N spacers) [3J, and will be 
disregarded in our study. One may, furthermore, show 
that any (time-reversal symmetric) quadratic cross term 
involving m, Q, and <f) does not respect the symmetries 
of our device geometry (keeping in mind that <f> — > —(f) 
under time reversal and m is a pseudovector under im- 
proper rotations). In particular, an interaction of the 
type cos(<^+rm z ) [10] is forbidden in our geometry. Thus 
neglecting interactions of m, Q, and (j) beyond quadratic 
order, the free energy remains uncoupled: 

F[m, Q, 4>] = F[m] + F[Q] + F[<f>] , (3) 

where F[m] = VKM 2 J2 - VM • H a , F[Q] = -QV, and 
F[(j)] = — Ej cos (j). The sign of the anisotropy constant, 
K, defines an easy plane or easy axis and is determined by 
the geometry of the device and crystalline anisotropies. 
H a is an applied external magnetic field. 

The LLG equation of motion of the magnet is now 
complemented with interactions that are quasistationary 
(i.e., first order in frequency), up to quadratic order in the 
components of m, preserving the magnitude of m, and 
consistent with the structural symmetry of the device: 

m = — 7m x H + am x m 

+ {pQ + \(f>)m x z x m + [yQ + K0)m x z . (4) 

Hereafter, we are focusing on the spin-valve case, 
Fig. [ljb), where the phenomenological coupling coef- 
ficients /1, A, v 1 and k may be taken to be angle- 
independent constants and z = m s . [For the SOI device, 
Fig.[T](a), structural symmetries dictate these coefficients 
to be odd functions in m z .] Constants \x and v charac- 
terize the strength of the coupling between the magnet 
and the total current Q. Similarly, the strength of the 
coupling between the magnet and the dynamics of the 
superfluid condensate <f> is characterized by A and k. To 
the reader familiar with spin valves [3J, Eq. Q is remi- 
niscent of the Landau-Lifshitz-Gilbert equation with the 
so-called Slonczewski and field-like torques, respectively, 
added on the second line of the right-hand side. In this 
case, sketched in Fig. [ijb) , current is spin polarized by 
passing through the fixed magnetic layer. The resulting 
spin-polarized current impinging on a free ferromagnet 
induces torque due to conservation of angular momen- 
tum. In the case of a single magnetic layer with SOI, 
Fig. [ija), a spin torque is generated via SOI inside this 
layer itself [15]. Because the leads in our system are 
superconducting, we additionally generate a torque as 
a result of the dynamics of the superfluid condensate. 
Loosely speaking, the torque induced by both currents, 
normal current and supercurrent, through the junction 
produce two channels for driving magnetization dynam- 
ics (and thus two sets of terms, as compared to the usual 
normal- metal spin torques). Appropriately, above the 
critical temperature of the superconductor, we expect to 
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recover the normal-metal limit, in which torque is gener- 
ated by the ordinary Ohmic current alone. 

The reaction of the current and superconducting phase 
dynamics to the magnet are not captured by the Joseph- 
son relations, Eq. ([2]), which would not be consistent with 
Eq. Q . One must extend Eq. ^ to include the pumping 
terms satisfying Onsager reciprocity, in order to obtain 
equations of motion for our coupled system that obey 
microscopic time-reversal symmetry [TT]. Because the 
magnet flips under time reversal (upon invoking Onsager 
symmetry), one must additionally use the symmetries 
of the structure to relate the time-reversed state to the 
original. After straightforward manipulations, that are 
analogous to Ref. |16j for normal junctions, we construct 
the following equations in lieu of Eq. ^ : 

Q = — [E j sin 4> — S(\m • m x z + urn ■ z)l + -r— , 

n 2e 

2e 

<j) = — [V — S(pm • m x z + ism ■ z)] — pQ , (5) 

where S = VM s /j is the total spin angular momentum of 
the ferromagnetic layer. These equations of motion now 
include both normal and superfluid pumping, which are 
Onsager reciprocal to the driving effects introduced in the 
generalized LLG equation, Eq. Q. Our theory includes 
two types of pumping as a result of the non- Ohmic rela- 
tionship between current and voltage. The term with co- 
efficient p causes current to drag phase across the device; 
p is a measure of the viscosity between the current and 
superfluid condensate. Although this term is not needed 
for consistency with Onsager reciprocity, we will see that 
it would generally have to be included in order to satisfy 
the second law of thermodynamics. We could also im- 
mediately notice that the coefficients p, is, and p should 
vanish in the limit of large superconducting reservoirs, 
recovering the ordinary ac Josephson effect (as expected 
based on the gauge invariance) . Keeping these terms, on 
the other hand, would capture finite-size (mesoscopic) 
properties of the superconducting layers, which are of 
secondary interest to our ends. 

We may write the equations of motion in a dimension- 
less form by measuring time, magnetic field, charge, volt- 
age, and conductance in units of S/Ej, Ej/jS, 2eS/h, 
Ejh/2eS, and S(2e/ti) 2 , respectively: 

m = m x H + am x m + </>(Am x z x m + Km x z) 

+ Q(pm x z x m + ism x z) , 
Q = sin (f> — Am • m x z — Km ■ z + a(j) , 
<f> =V — pm • m x z — vm ■ z — pQ . (6) 

Additionally, allow us to absorb a factor of VAf s 2 / Ej into 
the anisotropy constant, such that the free energy for the 
magnet reads F[m, Q, 4>] = Ej(Km 2 z /2 — m H a - QV — 
cos <p) . Under time reversal, the terms with coefficients 
v and A reverse sign in the LLG equation. Because rh 



does not change sign, these are dissipative. Likewise, 
the terms with coefficients p and k do not reverse sign 
and are thus nondissipative. a is a dissipative coefficient, 
therefore p is as well. 

Let us try to understand the microscopic origin of the 
dissipative terms in our theory. Consider momentarily 
only the RSJ subsystem: when enough energy is sup- 
plied (either thermally or by a bias), quasiparticles are 
able to overcome the superconducting gap and transport 
through the junction. Normal scattering of quasiparti- 
cles across the junction causes Ohmic resistance. Like- 
wise, consider an isolated precessing ferromagnet. This 
is microscopically described by a coherent magnon state 
that can decay into phonons and incoherent magnons, 
processes which macroscopically give Gilbert damping. 
In the case of a metallic ferromagnet, the additional de- 
cay channel into the electron-hole continuum enhances 
further its Gilbert damping. Upon coupling these sub- 
systems, energy is shared by the entire structure. Like- 
wise, dissipation by microscopic mechanisms underlying 
Ohmic conductance and Gilbert damping can give rise 
to a dissipative (viscous) energy transfer between ferro- 
magnetic and superconducting layers, as parametrized by 
new dissipative coefficients is and A. Phenomenologically, 
therefore, we may expect a , a, and p to bound is and A, 
which is indeed verified below. 

In the RSJ model, Eq. if cjy is static, we are in 
a superconducting (S) state because only dissipationless 
current is passing through the junction. Likewise if 4> is 
not constant, the circuit must have a finite voltage drop. 
This is called a resistive (R) state. Notice that in our 
generalized model, Eqs. ([6]), a choice of dynamics that 
leave (f> static can still generate dissipative current due to 
magnetic pumping. We will, nonetheless, keep refering to 
the static and dynamic states of </> as the superconducting 
(S) and resistive (R) states, respectively, even though this 
terminology is, in general, abusive, in the presence of the 
new spin-torque/pumping terms in Eqs. 

We distinguish between two regimes governed by the 
superconducting coherence length £. When £ is smaller 
than the width of superconducting terminals, the bulk 
properties of the superconductors will be largely detached 
from physics at the interfaces. Thus for large supercon- 
ducting reservoirs, a change in phase difference cannot 
be induced by transport through the junction. We ex- 
pect the corresponding coefficients p, is, and p to scale 
inversely with the volume of the smaller of the super- 
conducting layers then; these are representative of meso- 
scopic effects, as has already been inferred above. Be- 
cause charge is a conserved hydrodynamic quantity, on 
the other hand, there is a length at which the dynamics 
at the interface decouple from the current in the bulk. In 
particular, k, A, and a should not depend on the size of a 
superconducting reservoir; these coefficients parametrize 
the properties of the Josephson junction itself and are 
thus of central interest to us. Dynamic properties of 
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mesoscopic junctions with nonzero v and p terms in 
Eqs. ([6| are discussed in the Supplementary Text, where, 
in particular, we point out a reentrant behavior for the 
resistive (R) state as a function of the applied current I. 

For our analysis, it is convenient to consider dc current 
biasing (setting p, v, p to zero), Q = I: 

m = m x H + am x m + </>(Am x z x m + nm x z) , 
(Tcb = I — sin tf> + Am • m x z + Km • z . (7) 

In what follows, we treat a system where the applied 
magnetic field is along the axis of symmetry, H = H a z, 
and K is positive (which is generically the case for 
films with magnetostatic energy dominating over crys- 
talline anisotropy). Thermodynamic self consistency 
of our theory requires for the dissipation power P = 
(Ej/S)(ath 2 — 2Arh-mxz^+(Tf ) > 0. This bounds our 
phenomenological constant A as A 2 < aa (while, clearly, 
a > and a > 0). To proceed with the nonlinear stabil- 
ity analysis, it is important to notice that, according to 
Eqs. Q, the dynamics of m z and 4> decouple from the 
transverse magnetization (m x ,m y ), which can, in turn, 
be expressed in terms of (m z , <j>) [see Eq. ( 10 1 in the Sup- 
plementary Text]. 

When the current is below the critical current, I < 1, 
one can show that there are three stable fixed points: po, 
ao, and oq which correspond to a Josephson 0-junction 
(defining a junction with \(j>\ < ir/2 to be in the "0 phase" 
and 7r/2 < \<fr\ < tt in the u ir phase") and magnetic di- 
rection parallel, antiparallel, and away from the z axis, 
respectively. In all these states (j> is fixed by the applied 
current such that sin <j> = I. As indicated by our stabil- 
ity diagram, Fig. [2j the state of our device is determined 
by the applied magnetic field. When \h a \ < 1, where 
h a = H a /K, m z — h a and m is fixed by initial condi- 
tions. By applying a sufficiently large external magnetic 
field, \h a \ > 1, oq is annihilated under a saddle-node bi- 
furcation [17] , and the sole stable state is po or ao for pos- 
itive or negative applied field, respectively, pinning the 
magnet along the z axis. A full linear stability analysis 
is discussed in the Supplementary Text where we note, 
specifically, that the dissipation power bound precludes 
the existence of a 7r-junction. 

If the current exceeds its critical value, I > 1, the su- 
perconducting phase and z component of the magnet be- 
come dynamic. This disappearance of all the fixed points 
is an infinite-period bifurcation [17 . Because no fixed 
points exist, the Poincare-Bendixson theorem implies any 
closed orbit on the cylinder, parameterized by m z and cj>, 
is periodic and must go around the circumference of this 
cylinder. Supposing the frequency of this periodic motion 
is fij, m z may be written as a constant plus terms pe- 
riodic in f2j. Likewise, we may express <f> = nfljt (with 
nonzero n £ Z) plus terms periodic in Slj. Therefore, 
the characteristic frequency of the system is given by the 
time average of <fi. Upon solving the equation of mo- 
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FIG. 2. Stability diagram as a function of the current, 7, 
and applied magnetic field h a . A = —0.1, p, v, k = 0, K = 1, 
a = 1, and a = 0.1. h labels the Hopf bifurcation (solid lines), 
i labels the infinite-period bifurcation (long-dashed lines), and 
d labels the saddle-node bifurcation (short-dashed lines). 



tion for transverse component of the magnet [Eq. ( 10 ) in 
the Supplementary Text], we find it undergoes rotations 
at frequency n(X/a)VLj that are superimposed with VLj 
oscillations. Therefore the magnet in general undergoes 
quasiperiodic motion, a state we label q. 

To determine the full expression for flj when I > 1 
would require solving the system of differential equations 
([7]). For simplicity, consider the limit of small A and k, 
so that we can neglect quadratic terms in A and k in 
Eqs. ([7]). In this case, the characteristic frequency of the 
Josephson junction is given by the usual RSJ frequency 
ojj = V I 2 — 1/a [14] . In region q of our stability di- 
agram, Fig. [2] m z oscillates with frequency uij around 
the average value sign(J)(A/a — k)ujj + h a /K. Near the 
point |sign(/)(A/a — k)ujj + h a /K\ = 1, a Hopf bifurca- 
tion [T7] (labeled h) is induced wherein the quasiperiodic 
orbit disappears and the magnet is parallel or antiparal- 
lel to the z axis, labeled pu and apt. respectively, and the 
phase is dynamic. We anticipate the higher-order cou- 
pling in A and k to modify the frequency dependence on 
current. Furthermore, we expect that, near the line defin- 
ing the Hopf bifurcation, there exists a phase of bimodal 
stability wherein the magnet can orient along the z axis 
or precess quasiperiodically, subject to the initial condi- 
tions. This is a natural consequence of the reciprocity of 
current-driven magnetic dynamics and pumping and per- 
sists even in the absence of any superconductivity (i.e., 
Ej = 0). Details of these rich coupled nonlinear dynam- 
ics are, however, beyond the scope of the present paper. 

In summary, we have introduced a model of S|F|S and 
S|F|N|F|S heterostructures coupling the dynamics of the 
magnets to that of the superconductor via a Rashba SOI 
in single-layer junctions and via magnetic misalignment 
in spin-valve junctions. We expect such structures to be 
highly adaptable to uses in spintronics due to the versa- 
tility with which one can in principle influence both the 
magnet and superconductor. Chaos in ferrites and mag- 
netic thin films is often attributed to spatially nonuni- 
form magnetizations ]3] . Perhaps a simpler route towards 
chaos in our model is by applying a magnetic field per- 



pendicular to the axis of cylindrical symmetry. As a re- 
sult, the dynamic equations become three dimensional 
and thus no longer restricted by the Poincare-Bendixson 
theorem to periodic orbits or fixed points. 
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Supplementary Text 



Decoupled Junction. — In the special case where mesoscopic effects dominate (//, v, p ^ but A, k—0) in Eqs. (6), the 
current-biased magnetic and superconducting dynamics decouple. We take this opportunity to recall the properties 
of magnetic spin valves and the RSJ model of superconductors, to which the decoupled equations map. Ignoring A 
and k, Eqs. (6) become simply 

rh = m x H + am x rh + I(fim x z x m + vm x z) , cr<f> = I — sin (f> . (8) 

The equation of motion for the magnet is thus the LLG equation for a spin valve, including Slonczewski (/i) and 
field-like (y) torques, in the case that a fixed magnetic layer points along the z axis. The superconductor is described 
by the RSJ model with zero capacitance. Appropriately, we find that the dissipation power depends only on the 
dissipative constants: P = (E 2 /S){am 2 + 2vm ■ zQ + a<j) 2 + pQ 2 ) > 0. Because a > 0, a > 0, and p > 0, our 
phenomenological constant is bounded, v 2 < ap, as anticipated. There are three possible stable states of the 
current-biased magnet in the presence of a static field in the z direction: pinned parallel to the z axis, antiparallel 
to the z axis, or precessing around the z axis, labeled p, a, and o, respectively. A pinned state is stable when 
\(fi/a — v)ljK + h a \ > 1. If \(fi/a — v)I/K + h a \ < 1, the magnet precesses at frequency u>m = pl/a. The 
corresponding stability diagram with Hopf bifurcation lines is shown in Fig. [3j In the dimensionless form of the RSJ 
description, when — 1 < I < 1, the junction is in the S state and the phase is fixed at <f> = sin -1 /. When the current 
is raised beyond the critical current, I > 1, the Josephson junction is in the R state and <fr oscillates with frequency 
ujj = \/ 1 2 — l/o~. For the RSJ model, a n junction is trivially impossible: \(f>\ cannot access values between tt/2 and 
7T. The inset of Fig. [3] displays the well-known phase diagram of the RSJ junction. 
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FIG. 3. Stability diagram as a function of the current and applied magnetic field of the decoupled magnet, fj, — —1.5, v, A, k = 0, 
and K — 1. p and a label the parallel and antiparallel states of the magnet, respectively. Inset: decoupled Josephson junction. 
The S state (unshaded) and R state (shaded) are separated by the line 1 = 1. Solid line is the value of <p for a junction and 
dashed for the unstable -n junction. 



Stability Analysis. — We enumerate all of the fixed points of the junction discussed in the Main Text and calculate 
their associated stability. The equations of motion for m z and <fi decouple from the transverse dynamics, according 
to Eq. (7): 

I -sine/) - X(H a - Km z )(\ - m|) , g . 
a — (Xk + kA)(1 — m 2 ) 

where A = (A — an) /(! + a 2 ), R= (k + aA)/(l + a 2 ), and a. = a/(l + a 2 ). The solution for the transverse components 



m z = (1 — m 2 ) a(H a — Km z ) + X<j) 
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of the magnet, to = m x + im y , as a function of m z and (f> are given, in turn, by 



exp 



Ac/>- 



1 1 - m 2 
- In 

2 1 + to. 



(10) 



with (/? determined by initial conditions. Consequently, the fixed points of the equations of motion for m z and <fi, 
Eq. ([9|, immediately determine the state of the full system. When the current is below the critical current, 7 < 1, 
there are at most six fixed points for the (m z ,(f>) dynamics [Eqs. Q]: 



(m z ,(j)) = {(±l,sin 1 1) , (±l,7r-sin 1 1) , (h a ,sm 1 1) , (h a ,n 



(11) 



We henceforth label the first four fixed points as po, ao, p-x and a„. At these points, the magnet points parallel (p) or 
antiparallel (a) to the z axis, and the Josephson junction is a or it junction. We refer to the last two fixed points 
in Eq. (Ill (which are only possible when \h a \ < 1) as oq and o n , with m z determined by the ratio of the applied 
field to the magnetic anisotropy. When |ft Q | > 1, oq and are annihilated under a saddle-node bifurcation, and 
only the fixed points pinned along the z-axis, po, ao, p v and a ff , remain. At these fixed points, the superconducting 
phase and magnetic order decouple. The stability analysis for <j> then reduces to the RSJ model, resulting, therefore, 
in a junction. Thus, the only stable points, when \h a \ > 1 and I < 1, are po and a subject to the direction of the 
applied field. For intermediate values of the applied magnetic field, \h a \ < 1, po, ao, p^ and a w are all unstable fixed 
points and the system is at oo or o„. We can see that if a < (A« + kX)(1 — m 2 z ) [such that the denominator in Eq. ^ 
for <j) can become negative], the magnet is capable of sustaining the system in a state. Achieving a 7r junction 
for an optimal choice of K requires A 2 > aa/(l — to 2 ), in direct contradiction with the aforementioned dissipation 
power bound. Therefore, in this model, a ir junction is forbidden and, specifically when I < 1, the stable states of 
our system are po, ao, or oq, according to the value of h a . 



General Junction. — In the following, we analyze the properties of the general junction wherein we do not restrict 
any phenomenological parameters in Eqs. (6) to be zero. As previously, the transverse magnetization to = m x + im y 
decouples from the (m z ,(f>) dynamics: 



rh z = (1 — m 2 ) a(H a — Km z ) + \<f> + fil 



I - sin (f> - [X(H a - Km z ) - I(p,K + PA)] (1 - to 2 ) 
a — (\R + kX)(1 — m 2 ) 



(12) 



where fx = ([i — av) / {\ + a 2 ), v = (v + ct[i)/{l + a 2 ) and A = (A — an) /(1 + a 2 ), k = (n + aX) / (\ + a 2 ), a = a/(l + a 2 ), 
as before. One can show that the general solution for transverse components is (up to an arbitrary phase shift ip) 



exp 



1, 1 — m z 
lilt + X<t> + - In — ^ 
2 1 + rn z 



(13) 



The fixed points in the (m z , (f>) plane are 

(m z ,4>) = {(i^sin^ 1 /) , (±1, 7T - sin" 1 /) , ((fi/a - v)I/K + h a , sin" 1 7') , ((/x/a - v)I/K + h a , it - sin" 1 1') } , 
where we have introduced 



(14) 



/' = I [1 + (fi/a)X (1 - to 2 )] (15) 

with rh z = (fi/a — v)I / K + h a that itself depends on the current bias I. At the first four fixed points, the magnet 
is pinned parallel or antiparallel to the z axis and can be either a or n junction. Hence, maintaining consistent 
language between the coupled and general junctions, we label these fixed points po, ao, p^, and a^. The final two 
fixed points [which are possible when |(/i/ 1 a — v)I j K + h a \ < 1] are labeled by oo and o w . These oo and o„ points are 
stationary in the (to z , cf>) plane but the transverse components of the magnet follow a circular orbit of radius yjl — m 2 
at frequency lum = [ilja. 

The salient differences between these fixed points and those found studying the fixed points of Eq. ([9| are in the 
properties of oq and o w . First, the transverse component of the ferromagnet is dynamic when \x ^ 0. Second, the 
static value of sin cf> is a nonlinear function of the current. This results in a change in shape of the boundary separating 
the S and R states of the superconductor: See, for example, Fig. |4j where the phase diagram develops a "foldover 
region." Consider the current increase at fixed magnetic field along the dashed line in Fig. |4j The system undergoes 
changes from (1) S to R, (2) R to S, and (3) S to R again. Unlike in a conventional Josephson junction, our model 
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FIG. 4. Separation of the S state (white) and R state (grey) of the superconductor by a nonlinear function denned by /' = 1. 
The parameters of this system are /i = — 1, A = 0.6, u, k = 0, K = a = 1, and a — 2. The 1, 2, 3 labels along the dashed line 
show the three places where the Josephson junction switches between superconducting and resistive states. 



has multiple values of the current for which the junction changes between superconducting and resistive states. Thus 
the junction has three 'critical currents.' Likewise at a particular fixed value of current, we can induce a change from 
R to S then S to R by increasing or decreasing the applied magnetic field. This has no analogy in the RSJ model. As 
a function of the applied current, a rich variety of the coupled dynamics generally emerges, as seen in Fig. [5j where 
we have plotted the stereographic projection of the magnetic direction. A detailed analysis of this motion will be 
addressed in future work. 




FIG. 5. Stereographic projection of the magnetization undergoing irreversible dynamics at different currents. Here, fi = O.f , 
A = 0.5, v, k = 0, K = 1, a = 1, and <j = 1; initially positioned at m x = f . Note that the scale is different between frames. 



